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ABSTRACT 

We have obtained a time series of 81 high-cadence circular polarization observations of the rapidly oscillating Ap star HD 24712 
with the new ESPaDOnS spectropolarimeter at CFHT. We used the high-S/N, high-resolution Stokes / and V spectra to investigate 
possible variation of the mean longitudinal field over the pulsation cycle in this roAp star. Our multiline magnetic field and radial 
velocity measurements utilized 143 spectral lines of rare-earth elements, attaining precision better than 13 G and 19 m s , respectively. 
A multiperiodic radial velocity variation with an amplitude of 40-136 ms~' is clearly detected at the known pulsation frequencies of 
HD 24712. At the same time, no evidence for pulsational changes of the magnetic field can be found. We derive a 3cr upper limit of 
10 G, or about 1% of the mean longitudinal field strength, for magnetic field oscillations in the upper atmosphere of HD 24712. The 
absence of detectable pulsational variability of the magnetic field provides a valuable constraint for the interaction between pulsations 
and magnetic field in roAp stars and is compatible with the recent predictions of detailed theoretical models of stellar magnetoacoustic 
oscillations. 
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1. Introduction 

Rapidly oscillating magnetic Ap (roAp) stars are cool magnetic 
chemically peculiar stars exhibiting a remarkable combination 
of unusual surface properties. These objects have strong, organ- 
ised magnetic fields, most likely of fossil nature, and show pho- 
tospheric chemistry deviating far from that of the Sun and other 
cool, non-magnetic A-type stars. They show non-uniform dis- 
tributions of chemical elements, both laterally across their sur- 
faces, as well as vertically with height in their atmospheres. Most 
remarkably, they exhibit high-overtone, low-degree, non-radial 
p-mode pulsations with periods ranging from 6 to 21 min. 

Their strong magnetic fields play a central role in the mech- 
anisms responsible for pulsation mode excitation and selection 
(Balmforth 2001] Saio[2005 ) and in shaping the frequency spec- 
tra of roAp stars (Cunha 120061) . Pulsational perturbations ob- 
served in roAp stars occur in the magnetically-controlled outer 
stellar layers, which is why pulsations are aligned with the 
oblique quasi-dipolar magnetic field (Kochukhov 2004; Saio 
120051 ) instead of the stellar rotation axis as in all other types of 
non-radially pulsating stars. 

Recent time-resolved spectroscopic observations of pulsa- 
tions in roAp stars (Kochukhov & Ryabchikova 2001a, 2001b; 
Mkrtichian et al. I2U03I Elkin et al. [20051 Kochukhov [20061 
demonstrate that, in addition to the remarkable horizontal ge- 



ometry, the magnetoacoustic p-modes show outstanding vertical 
structure due to the combination of the rapid increase of pul- 
sation amplitude with height in the atmosphere and the short 
vertical wavelength of the pulsational fluctuations. The gen- 
eral picture of the atmospheric pulsational behaviour of roAp 
stars is determined by the propagation of magnetoacoustic waves 
through the distinct layers of the chemically stratified atmo- 
sphere (Ryabchikova et al. |2002). Pulsations are weak or non- 
detectable in the lower atmosphere probed by lines of light 
and iron-peak elements. On the other hand, rare-earth elements 
(REE) are concentrated in clouds located high above the pho- 
tosphere, in the layers where pulsation amplitudes reach up to 
several kms -1 . 

An important new insight into the complex physics of pulsa- 
tions in roAp stars can be achieved by the analysis of interaction 
between the p-modes and the magnetic field. The first search for 
variations of the mean longitudinal magnetic field, (B z ), over the 
pulsation cycle was conducted by Hubrig et al. (2004 ). These au- 
thors used low-resolution time-resolved circular polarization ob- 
servations of six roAp stars, and measured (B z ) using hydrogen 
lines and unresolved blends of metal lines. Hubrig et al. (2004) 
failed to detect magnetic field variation above the noise level of 
their analysis, which was 40-100 G. 
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Kochukhov et al. (2004b) and Savanov et al. ([2006) obtained 
high-precision measurements of the mean magnetic field mod- 
ulus over the 12-min pulsation period of the bright roAp star 
y Equ. The resolved, Zeeman-split profile of the Fe n 6149.25 A 
line was used in both studies. Neither of the authors detected 
any pulsational fluctuation of the field strength, deriving an up- 
per limit of 10 G for the possible field variation. 



The low-resolution spectropolarimetry and the field modu- 
lus measurements using Fen lines have limited diagnostic ca- 
pabilities because both types of magnetic field monitoring tech- 
niques probe low atmospheric layers where pulsation amplitudes 
are small. Taking guidance from the highly successful radial ve- 
locity analyses, an investigation of the magnetic field oscilla- 
tions in roAp stars should focus on individual REE lines which 
show conspicuous pulsational variations. The first study of this 
kind has been carried out by Leone & Kurtz (2003). They ob- 
served several Nd m lines in y Equ using a circular polarization 
analyzer and a high-resolution spectrometer. Based on only 18 
time-resolved spectra, Leone & Kurtz (2003 ) announced the dis- 
covery of the pulsational variation of (B z ) with amplitudes 1 10- 
240 G although curiously with discrepant phases of magnetic 
maximum for different Ndm lines. Subsequently, Kochukhov 
et al. d2004al l obtained independent high-resolution polarimetric 
observations of y Equ. Their analysis combined magnetic signa- 
tures of 13 Ndm lines and was based on more than 200 spectra 
collected over three nights. Kochukhov et al. (12004a) were un- 
able to detect longitudinal field variability over the pulsation cy- 
cle in y Equ exceeding 40 G. They suggested that the field vari- 
ation claimed by Leone & Kurtz (2003 ) is spurious and likely 
results from the neglect of blending of the few REE spectral fea- 
tures employed in their magnetic measurements. 

The contradictory results obtained in the previous searches 
for rapid magnetic field oscillations in y Equ call for an exten- 
sion of time-resolved spectropolarimetric studies of REE lines 
to other roAp stars. We have therefore obtained a new set of 
time-resolved, wide wavelength coverage Stokes / and V spectra 
for the well-known roAp star HD 24712 (HR 1217, HIP 18339, 
DO Eri). 

HD24712 was identified as a roAp star by Kurtz d!982| l. 
The recent extensive Whole Earth Telescope photometric cam- 
paign revealed eight pulsation modes with frequencies in the 
range of 2.6-2.8 mHz (Kurtz et al. 2005). Ryabchikova et 
al. ([1997) performed detailed model atmosphere and chemical 
abundance analysis of HD 24712. They inferred T e g = 7250 K 
and log g =4.3, and found that the atmosphere of this roAp star 
is rich in REEs but relatively poor in the iron-peak elements. In 
their study of the evolutionary state of magnetic chemically pe- 
culiar stars Kochukhov & Bagnulo (2006) derived the following 
fundamental characteristics of HD 24712: mass M- 1.55 M Q , 
luminosity L-7.4 L G and an age « 10 9 yr, corresponding to 
about 50% of the main sequence lifetime. Bagnulo et al. ( 1995) 
demonstrated that a dipolar magnetic field topology with the po- 
lar strength B p — 3.9 kG and the angle B = 150° between the ro- 
tation and magnetic axes provides a good fit to both the lon- 
gitudinal field curve of HD24712 and to the broad-band lin- 
ear polarization measurements. Combining magnetic field and 
line strength observations of HD 247 12 obtained over more than 
three decades, Ryabchikova et al. (2005) determined a rotation 
period of 12.45877 d. 

2. Observations and data reduction 

Eighty-one new Stokes V spectra of HD24712 were obtained 
on 9 January 2006 UT using the ESPaDOnS spectropolarimeter 
at the Canada-France-Hawaii Telescope. Observations covered 
roughly 3 hours, starting at HJD = 24 53744 .70540. According 
to the ephemeris of Ryabchikova et al. (2005 ), our time-resolved 
monitoring of HD 24712 corresponds to the rotation phase 0.9, 
which is close to the positive extremum of the longitudinal mag- 
netic field variation and to the phase of maximum non-radial os- 
cillation amplitude. Observations were obtained in polarimetric 



mode, with a nominal resolving power A/AA = 65000, and cover 
the spectral region 3700-10500 A. The details of the ESPaDOnS 
are reported by Donati et al. (in preparation^ In particular, 
both results of the commissioning tests and our own velocity 
measurements using telluric absorption lines in the spectra of 
HD24712 indicate that the typical relative spectral stability of 
the instrument is better than 10-15 m s _1 for the frequency range 
corresponding to roAp pulsations. 

The 81 Stokes V spectra of HD 24712 were computed from 
162 individual circular polarization exposures. Between subse- 
quent exposures, the Fresnel rhomb quarter-wave retarder was 
rotated by 90° (i.e. alternately to position angles of +45° and 
-45°) in order to interchange the positions of the ordinary and 
extraordinary beams on the CCD. The data were reduced using 
the Libre-ESpRIT package; each subsequent pair of CCD frames 
was combined into a single Stokes V spectrum with typical peak 
S/N of 1 15. Examination of the diagnostic null spectrum (Donati 
et al. 119971 ) indicates that the beam-switching technique has re- 
moved all spurious contributions to the Stokes V spectrum to 
below the noise level. We emphasize that this is the first time 
that this high-precision spectropolarimetric observing method 
has been applied to monitor rapid magnetic field variations of 
a roAp star. 

Individual circularly polarized exposures were of 20 seconds 
duration, and the readout using the CCD's fast readout mode re- 
quired 25 seconds. Taking into account time required for repo- 
sitioning of the waveplate between exposures, the average time 
required for acquisition of a single Stokes V spectrum was typ- 
ically 1.25 minutes. This effective time resolution of the time 
series is significant in comparison with the P = 6.0-6.4 min pul- 
sational variation of HD 247 12. Our numerical modelling shows 
that the pulsational amplitude inferred from such a time series 
is underestimated by 12% relative to observations with a much 
higher time resolution. The treatment of this effect must be in- 
cluded in the time-series analysis. 

3. Time-series analysis of radial velocity and 
magnetic field 



Ryabchikova et al. ( 2007]) have recently presented a detailed 
study of the spectroscopic pulsational behaviour of HD24712 
using high-resolution echelle spectra acquired at several large 
telescopes. This analysis demonstrated that, similar to other 
roAp stars, the pulsational radial velocity (RV) variation can 
be clearly detected only for REE lines and, marginally, for the 
cores of a few very strong lines of light elements (hydrogen 
Balmer lines, Can 3933 A). Taking these results into account, 
we have investigated the pulsational variation of the Stokes / 
and V profiles for REE lines. An initial list of the central wave- 
lengths and Zeeman splitting parameters was compiled for more 
than 200 lines based on the information from Ryabchikova et al. 
(120061 120071 and the VALD (Kupka et al. [T9991 and DREAM 
(Biemont et al. 1999 ) databases. This list is dominated by Prm, 
Nd ii and Nd in lines, but also includes Y n, La n, Ce n, Pr n, Sm n, 
Eun, Gdn, Tbin, Dyn, Dym, and Erm spectral features. The 
Stokes / and V profile of each line was examined to remove lines 
affected by blending. Radial velocities were determined for the 
remaining lines with the center-of-gravity method: 
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Fig. 1. Pulsational variation of radial velocity (upper panels) and longitudinal field (lower panels) for the REE lines in HD 24712. 
The left panels show the amplitude spectra of RV and (B z ) (thin line) as well as Fourier transforms of the residual variation (thick 
line) obtained by subtracting least-squares sinusoidal fit from the original time series. The vertical bars show pulsation frequencies 
of HD 24712 determined in the photometric study by Kurtz et al. (2005 ). The right panels illustrate velocity and magnetic field 
variation in the time domain. Individual REE and 5(B Z ) measurements (symbols) are compared with the best fit three-frequency 
model (thick line), which takes into account data sampling and phase smearing (see text). The thin line shows the same oscillation 
model neglecting these effects. 



where Aq is the laboratory wavelength, / and I c correspond to the 
line and continuum intensity, respectively, and c is the speed of 
light. 

The mean longitudinal magnetic field was estimated by inte- 
grating the Stokes V profiles of REE lines: 



(B z ) = - 



4nm e c J (v - (v»Vdv 
e Aoz f(h-I)dv' 



(2) 



where z is the effective Lande factor, m e is the electron mass, e is 
the electron charge, and v is the velocity coordinate. The errors 
of (v) and (B z ) were derived applying the standard error propa- 
gation rules to Eqs. (Q~|) and (|2]), and using the formal uncertain- 
ties of the Stokes / and V spectra provided by the Libre-ESpRIT 
spectral extraction package. 

Preliminary time-series analysis was carried out for each 
measured transition to exclude lines showing no RV oscilla- 
tion signal above 2cr in the 2.4-3.0 mHz frequency interval to 
which all known pulsation frequencies of HD 247 1 2 are confined 
(Kurtz et al. 120051 1. Four Tb in were also excluded from consid- 
eration because their pulsational radial velocity curves show a 
large phase shift with respect to the variation of other REE ions 
(Ryabchikova et al. [2007l >. This left us with 139 REE lines. Due 
to overlap of the blue echelle orders in the ESPaDOnS spectra, 
many lines could be measured twice, increasing the total number 
of spectral features included in the final sample to 168. 

We have constructed a high-precision measure of the pul- 
sational changes of RV and magnetic field by calculating the 



weighted average of the deviation of (v) and (B z ) from their 
mean values: 



(ai + bit)] /o-j(S) 



Zf =1 l/o*(S) 



(3) 



where N = 168, S corresponds to (v) or (B z ), o-j(S) are the 
respective errors and a,, b, represent coefficients of a linear least- 
squares fit to the full timeseries describing the mean and linear 
drift of each quantity. 

The amplitude spectrum and time dependence of the result- 
ing average 6(v) and 6(B Z ) are illustrated in Fig.[T] We have un- 
ambiguously detected pulsational variation of RV at the level of 
« 150 ms . Pronounced modulation of the RV amplitude due 
to beating of the oscillations with P « 6.05 and 6.30 min is evi- 
dent in our data. The first period coincides with the frequency V5 
(2755 yuHz, 6.049 min) of Kurtz et al. (12005 ), whereas the second 
one lies between v\ (2620 yuHz, 6.362 min) and V2 (2653 /uHz, 
6.282 min). Given that the frequency spectrum of HD 24712 is 
well known and remains constant over many years, we assume 
that the p-mode oscillations detected in our CFHT spectra are 
caused by the superposition of the three pulsation periods corre- 
sponding to vi, V2 and V5 of Kurtz et al. ([2005 )■ We keep these 
periods fixed in the following time-series analysis. 

Figure Q] shows that there is no evidence of the (B z ) varia- 
tion with any of the known pulsation frequencies of HD 24712, 
nor any other frequency up to the Nyquist limit of our data. The 
highest 6(B Z ) peaks reach 6-8 G, but none is statistically signif- 
icant. 



Table 1. Results of the multifrequency least-squares analysis of the variation of radial velocity and longitudinal magnetic field for 
REE lines in HD 24712. For each group of lines we list the number of measured spectral features, followed by the number of unique 
lines. Amplitudes are given in m s _1 and G for (v) and (B z ), respectively, while phases are measured in fractions of pulsation period. 



p 




Radial velocity 


Magnetic field 






(min) 


A (m s" 1 ) 


<p A (G) 




w 

T 












168(139) REE lines 












6.049 


1 38 7 


+ 3 Q 

x j.y 


0.850 ± 0.005 4.9 ± 2.4 


0.857 


± 


0. 


076 


6.282 


1 13.9 


+ 6.2 


0.861 ±0.009 6.6 


±3.8 


0.143 


± 


0. 


098 


6.362 


40 


+ 63 

X KJ.J 


0.120 + 0.024 5.1 ±3.9 


0.944 


± 


0. 


118 








56(47) Nd ii lines 












6.049 


1 80 


+ 6.0 


936 + o 005 5 2 


± 5.6 


0.399 


± 


0. 


.174 


6.282 


132.0 


+ 9.5 


0.971 ±0.012 19.9 


± 9.0 


0.317 


± 


0. 


072 


6.362 


41 3 


+ Q 7 


HQ + 037 8 


±9.1 


0.808 


± 


0. 


184 








44(34) Nd m lines 












6.049 




+ 4 3 


0.842 ± 0.005 4.8 


±3.2 


0.944 


± 


0. 


106 


6.282 


138.6 


+ 6.9 


0.851 ±0.008 7.4 


±5.1 


0.196 


± 


0. 


108 


6.362 


4^ 7 


-+- 7 1 

X / . 1 


0.148 ±0.026 10.3 


±5.2 


0.037 


± 


0. 


079 








17(16) Prm lines 












6.049 


179.5 d 


: 8.6 


0.662 ± 0.007 16.5 ± 


7.6 


0.859 


± 


0. 


075 


6.282 


173.7 =t 


: 13.5 


0.663 ±0.012 18.4 ± 


12.1 


0.902 


± 


0. 


105 


6.362 


73.8 d 


: 13.8 


0.975 ± 0.030 8.7 ± 


12.2 


0.754 


± 


0. 


227 








15(12) Smn lines 












6.049 


213.5 d 


: 10.5 


0.923 ± 0.008 10.5 ± 


9.2 


0.829 


± 


0. 


143 


6.282 


191.4d 


: 16.9 


0.966 ± 0.014 30.0 ± 


14.9 


0.678 


± 


0. 


078 


6.362 


26.3 d 


: 17.2 


0.419 ±0.104 18.4 ± 


15.1 


0.474 


± 


0. 


130 



HD 24712 exhibits no pulsational variation of the disk-averaged 
line of sight magnetic field component. The 3cr upper limit for 
pulsational field changes is 7-12 G, depending on the pulsation 
mode considered. 

We have examined separately the average RV and (B z ) for 
several groups of REE ions. Table[TJpresents results of this anal- 
ysis for the four groups of REE species with > 10 unique lines 
each (Prm, Ndn, Ndin, Smn). The time-series analysis is less 
precise for individual ions, but it generally confirms results ob- 
tained with all REE lines. In particular, none of the investigated 
REE species shows statistically significant pulsational changes 
of (B z ), although each of them shows RV pulsations at >20<x 
level. 

The REE line analysis was complemented by the study of 
(B z ) measured for the Least Squares Deconvolved Stokes / and 
V spectra (LSD profiles, see Donati et al. U9971 l. This diagnostic 
is based upon strong lines of light and iron-peak elements and, 
therefore, probes conditions in the lower part of HD 247 1 2 atmo- 
sphere. The amplitude spectrum of the LSD (B z ) measurements 
is presented in Fig. [2] No variation above 13-21 G (3cr limits) 
is present at the frequencies detected in the REE line radial ve- 
locities. The average longitudinal field obtained with the LSD 
technique is 1022 + 30 G, which is in good agreement with the 
known rotational modulation of (B z ) in HD 247 12 (Ryabchikova 
etal. l2005l[2007t . 

4. Discussion 

We have obtained the first time-resolved high-resolution circular 
polarization observations of the roAp star HD24712. Our aim 
was to search for the variation of the stellar magnetic field with 
the pulsation cycle. Unlike all previous time-resolved spectropo- 
larimetric studies of roAp stars (Hubrig et al. 2004; Leone & 
Kurtz [20031 Kochukhov et al. 2004a), we have employed the 
beam-switching technique to minimize spurious polarization. 
The resulting phase smearing was accounted for in the time- 



We have used a least-squares analysis in the time domain 
to determine RV amplitudes and to infer accurate upper limits 
for the magnetic field variation. For high-cadence data in which 
the time sampling is much shorter than the pulsation period, the 
expression 

*S=^AiCos — +2n<fA (4) 

!=1 ' ' ' 

(where A,- and tpi are the amplitude and phase of oscillation with 
period P,) is suitable to approximate the variability of pulsating 
stars. However, the phase smearing caused by the longer effec- 
tive time sampling used in this study is non-negligible. We have 
explicitly taken the data sampling effects into account by inte- 
grating and averaging Eq. <j4j for the time intervals correspond- 
ing to the individual polarimetric subexposures. Applying trivial 
trigonometric relations, one can obtain the following expression: 



— ^AA . (nAt\ ( nt^ + nt^ \ 

5S — > Sin COS h 27TIZ>, 

tf^Af \p t ) \ Pi Vl ) [5) 

xcos ^^ 45 ~ 7rr45 ), 



where f +45 and f~ 45 are mid times of the subexposures with dif- 
ferent retarder orientation and At — 20 s is the exposure time. 

Results obtained by fitting Eq. (0 to the (v) and {B z ) varia- 
tions are summarized in Table[TJand illustrated in Fig.Q] The RV 
pulsations have amplitudes 40-136 m s~' and are detected at the 
6-36<x level. The rms deviation between observations and the 
least-squares fit indicates an uncertainty of 19 m s _1 for a single 
multiline RV measurement. The corresponding precision of the 
multiline magnetic measurements is 13 G. Fig.[TJshows that the 
amplitude of the residual oscillations does not exceed 5-10 m s~' 
after prewhitening observations with the three-frequency model. 
At the same time, the formal (B z ) amplitudes are 5-6 G, detected 
at the 2<t confidence level or less. Therefore, we conclude that 




Frequency (mHz) 

Fig. 2. Amplitude spectra of the LSD (B z ) measurements for 
HD 24712. Fourier transform of the original (B z ) time series 
(thin line) is compared with the residual amplitude spectrum 
(thick line). The vertical bars show photometric pulsation fre- 
quencies of HD 24712. 



series analysis. To boost precision further, we have combined os- 
cillation signal in 139 REE spectral lines, thus significantly im- 
proving the sensitivity to the magnetic field variation in compar- 
ison with previous studies. Variability with the known pulsation 
frequencies was found for the RV obtained from averaged REE 
lines as well as for individual ions. However, no magnetic field 
variability above 7-12 G could be detected. Since our observa- 
tions were obtained at the rotation phase of magnetic maximum 
of HD 24712, when (B z ) » 1 kG, we infer the upper limit of the 
relative pulsational changes of the field strength: 6B/B ^ 10 2 . 

The null result of our search for the pulsational magnetic 
field variation in HD24712 agrees with the 2^1% 6B/B upper 
limit obtained by Kochukhov et al. (2004a) for y Equ. Our study 
thus strengthens the view that no significant field variability in 
roAp stars should be expected on pulsational timescales. 

Hubrig et al. (2004) published a preliminary theoretical es- 
timate of the expected pulsational modulation of (B z ) in roAp 
stars. These authors claimed that field variations at the level of 
1% to 14% for different roAp stars may be expected. For the 
pulsation period and the radial velocity amplitude observed in 
the upper atmosphere of HD 247 12, magnetic changes should be 
about 1% of (B z ) or 10 G. This level of the magnetic field varia- 
tions is marginally inconsistent with our conservative upper limit 
of 7 G derived for the highest-amplitude pulsation mode. At the 
same time, Hubrig et al. (2004) predict that the rapid magnetic 
variations of y Equ should reach 10% of the field strength, which 
is clearly ruled out by the observational results of Kochukhov et 
al. (I2004ab . 

Saio (2005 ) has pointed out that in their estimate of the mag- 
netic field variation Hubrig et al. (2004) have erroneously disre- 
garded the effect of the horizontal component of the pulsational 
displacement. Based on the detailed non-adiabatic theoretical 
model of the magnetoacoustic oscillations in roAp stars, Saio 
(2005 ) suggested that the amplitude of the magnetic field varia- 
tion should not exceed 6B/B~ 10~ 5 or ~0.01 G for a star with 
a kG-strength magnetic field. This theoretical prediction is con- 
firmed by our spectropolarimetric monitoring of HD 24712 and 
by the time-resolved magnetic observations of y Equ presented 
by Kochukhov et al. (20d4at. 
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